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Abstract 



In the context of nonlinear scattering, a continuous wave incident onto 
a nonlinear discrete molecular chain of coupled oscillators can be partially 
absorbed as a result of a 3-wave resonant interaction that couples two 
HF- waves of frequencies close to the edge of the Brillouin zone. Hence 
both nonlinearity and discreteness are necessary for generating this new 
P^ ' absorption process which manisfests itself by soliton generation in the 

medium. As a paradigm of this nonlinear absorption we consider here 
the Davydov model that describes exciton-phonon coupling in hydrogen 
bonded molecular chains. 

t^j- ' 1 Introduction 

The scattering of waves becomes extremely rich when nonlinearity comes into 
^O play, and one of its most fundamental effect is soliton generation which results 

in energy localization and propagation. For instance an incident radiation on a 
two-level medium at the resonant frequency can be totally transmitted, instead 
of being absorbed, a property described in 1, as self induced transparency. It 
results from the nonlinear coupling of radiation with medium population, a 
mechanism which generates the solitons, vectors of energy transmission. 

Other interesting processes of nonlinear wave scattering are the two-photon 
i ^ ' propagation, second harmonic generation and stimulated Raman scattering. In 

that last case, the nonlinear interaction induces (laser) pump depletion and 
Vh ' phase effects result in Raman spike generation (short duration pump repletion) 

0. Although the Raman spike generation is not a solitonic effect, here also the 
nonlinearity is the fundamental tool JHJ- 

Recently discovered, the nonlinear supratransmission 01 E] is a nonlinear 
scattering process where a totally reflective medium switches to high transmis- 
sivity above some threshold because of nonlinearity. The bifurcation takes its 
origin in a nonlinear instability 6 which has a wide field of application in generic 
nonlinear evolutions [Jj. 

Particularly interesting is the application of this concept to the scattering 
of a continuous-wave laser beam incident onto a Bragg mirror which switches 
from total reflection to transmission by means of gap soliton generation [5J 1!|] ■ 



It has been demonstrated that the switch is a manifestation of nonlinear supra- 
transmission which has in particular allowed to compute an analytic explicit 
expression of the intensity flux threshold |1U) . 

We are interested here in the concept of nonlinear absorption where a me- 
dium, transparent to incident radiation in the linear regime, can become ab- 
sorbant under contribution of nonlinearity which is due here to coupling of waves 
of different nature. Such is the case with the Davydov model :TT] describing 
the coupling of high frequency optical phonon waves to low frequency acoustic 
waves. A generic field of application of the Davydov model is the cristallin 
acetanilide where the acoustic wave represents the hydrogen bond between ad- 
jacent molecules, while the optical phonon wave represents the C=0 stretching 
(amide-I mode). This model has been widely studied in its stationnary limit 
where the hydrogen bond is supposed to have a time response much longer than 
the high frequency C=0 stretching. In that situation the resulting simplified 
model is the nonlinear Schrodinger equation that has been used to interprete 
the anomalous IR-absorption band by soliton generation |12| . 

The Davydov model has been recently reexamined in the context of reso- 
nant two-wave interaction where a rigorous multiscale analysis transforms the 
discrete model to a continuous integrable partial differential equation which 
keeps the fundamental wave coupling effect and possess multi-soliton solutions 
|13|. This two- wave resonant interaction process takes into account an inci- 
dent high-frequency optical phonon wave that resonates with the low-frequency 
acoustic wave as soon as the HF-wave group velocity equals the LF-wave phase 
velocity. This is the Benney criterion J3] that has revealed its efficiency through 
numerical simulations of the Davydov model |15| . 

The two- wave resonant process assumes no backward propagation and thus 
neglects the multiple reflections due to the periodicity of the cristal. We consider 
here a three-wave interaction process that couples the incident optical phonon 
wave both to the acoustic excitation and the reflected optical phonon wave. 
We shall discover that an incident and a reflected high-frequency waves can 
cooperate resonantly with a low-frequency acoustic wave thanks to the discrete 
nature of the basic model. 

We will demonstrate that this scattering process allows for absorption of 
incident radiation by a purely nonlinear effect that generates a three- wave soliton 
in the medium. The nonlinear effect is effectively an absorption as a part of the 
incident energy flux is transfered to the medium. 

2 Three- wave interaction in the Davydov model 

2.1 Basic model 

Our starting point is the Davydov model JTj for the eigenstate a n (i) of the 
amide-I excitation (the corresponding dynamical variable represents the C=0 
stretching) coupled to the dynamical variable j3 n that represents the longitudinal 



displacement along the hydrogen-bond chain. It reads 

iha n = [E + W + x(ftn+i - /?«)] a„ - J(a n+1 + a n -i), (1) 

Mf3 n = K(/3 n+ i - 2/3 n + /3„_i) + X (\a n \ 2 - K-i| 2 ), (2) 

where M is the mass of the peptide group, K the spring constant of the hydrogen 
bond, E the energy of the C=0 stretching and W the the total energy of 
the peptide group displacements. The constant % is the coupling parameter 
and J measures the energy of the dipole-dipole interaction of C=0 stretching 
oscillations. Overdot stands for derivation with respect to the physical time T. 
Upon defining the dimensionless time t — JT/h, and the adimensional quan- 
tities 

*„ - a n — %= exp[- T (Eo + W-2J)t], Q n = ^(0 n+ i ~ Pn) , (3) 
J\JM J <J 

the system becomes 

idt^n + (*n+l - 2*„ + *„-i) = Q„*„ , (4) 

d 2 t Q n - V 2 (Q n+1 - 2Q n + Q n _ x ) = |*„ +1 | 2 - 2|*„| 2 + |*„_i| 2 , (5) 

which constitutes our basic example of a discrete nonlinear coupled wave system 
of equations. Note that we are left with a single constant, the adimensional 
sound velocity 

V=-v p , Wp = y— , (6) 

where v p the phonon velocity (cells per second), and that the coupling parameter 
X has been absorbed in the amplitude of ip n {t). 

2.2 Multi-scale expansion 

Following the multiscale expansion method |16| we assume a representation of 
the solution {"^(i), Q n (t)} as formal series in terms of a small parameter e 
where space-time dependences occur at a sequence of slow scales. Since it is 
not necessary to push the series at arbitrary order, we write explicitely the first 
relevant terms only as 

* n (t) =eip i0) (n , Xl ,--- ]t ,t 1 ,---) + e 2 ip ( - 1 \n a ,x 1 ,--- ; to, t u ■ ■ ■ ) + • • • (7) 
Qn{t) = eq {0) (n Q ,xi,- ■■ ;*o,*i,- ••) + e 2 q {1) (n Q ,xi, ■ ■ ■ ;t ,ti, ■ ■ ■ ) H (8) 

and the difference-differential operators must be understood as 

V T t - V± + ed xi + ■ ■ • , dt -> d t0 + ed tl +■■■ . (9) 

Hereabove the difference operators are defined as 

V+tf„ = tf n+1 - tf n , V~* n = 9 n - *„_i, (10) 



such that the second order difference appearing in Q is 

V 2 *„ = V+V-* n = tf n+1 + *„_! - 2* B . (11) 

The system ll 1151) at first order gives the linear equations 

A>v (0) =o, £o = id to +v 2 no , 

W 0) =0, L = dl~V^l . (12) 

The next order e 2 eventually reads 

£o^ (1) = -id tl ^ - (V+ +V m )d x M 0) +q {0 ¥°\ 

W 1} = 2d to d t J 0) +^ 2 (V+ + V- )fl Xl gW + V£><°>| 2 . (13) 

The method then works as follows. Once selected an explicit solution of the 
linear system 11120 . we express that the evolution (|13[) of the first order cor- 
rection {tph , qn '} must not produce secular growth. This explicitely furnishes 
the evolution of the fundamental {</?„ , g„ } in the slow variables xi and t\. 
The choice of the linear solution in the variables noi^o determines the physical 
problem under study. 

2.3 Selection rules 

To describe a 3-wave interaction process involving incident and backscattered 
optical phonon waves, we are led to select in the linear system (|12f) the solution 

^ (0) = Aix^t^e 1 ^ 1 " ^ 1 ^ + B(x u t 1 )e t{ - k2n °- UJ2to) , (14) 

g (0) = g(xi, ti)e^ Kn "- at ^ + g(xi, t x )e-^ Kn °' ato) , (15) 

with wave numbers k\ > and fc 2 > such as to ensure an optical phonon wave 
as a superposition of an incomming wave of amplitude A and a reflected wave 
of amplitude B. The acoustic wave can propagate in both directions, thus K 
can be of either signs. We have explicitely written the variables {x\,t\) in the 
slowly varying amplitudes A, B, and g, but of course they depend on all higher 
order variables (but not on the first order ones no and to). The above expression 
is a solution of l|12|) for the following dispersion relations 

K 
uii = 2(l-cosfci) , uJ2 = 2(l-cosfc 2 ), Q = 2V|sm— |. (16) 

Note that using the relation J2J and the multiscale derivative laws, we 
demonstrate that the first order enveloppe /3o of the longitudinal displacement 
along the hydrogen-bond chain f3 n is related to the enveloppe g(x\, t\) of Q n by 

Po(zi,h) = -i- — r-7^:9(xi,ti) (17) 

2xsm(il) 



then the behaviour of the enveloppe g{x\ 1 t\) will automatically gives the be- 
haviour of (3a{xi,t\). 

The resonant wave interaction results from a selection rule for the wave 
parameters, obtained by examination of the nonlinear terms that occur in (|13|) . 
namely 

i g Qe i(-k-2 + K)n a e -i(u-2+ri)t _|_ gq e i(-k 2 -K)n e —i(ui2—Ci)to 

Such terms will combine to either components of ip(°> and resonate with corre- 
sponding factors in the left-hand-side of (|13fl . The evolution of the envelopes 
will be then obtained by setting to zero the coefficients of resonating terms. 

Since the physical context is the resonant interaction of two HF-waves (op- 
tical phonon) with a LF acoustic wave, it implies a small value of K (near the 
center of the Brillouin zone) and large values of fci and hz, that is to say close to, 
but less than, the value n in order to achieve incident and reflected HF-waves. 
We are then left with the following selection rules. First when the nonlinear 
terms Ag and Bg combine respectively with B and A we get 

fcl + fc2 = 2tt - K, K>0, ui-u 2 = -tt, (18) 

Another possibility is to combine instead Ag and Bg with B and A to obtain 

fcl + fc2 = 2tt + K, K<0, wi-W2=ft. (19) 

We shall study these two cases and demonstrate that only the first one gives an 
instability that generates solitons, it is displayed on fig-Q] 

It is worth remarking that such a scattering process (involving 2n) is allowed 
for by the presence of an exponential in the discrete variable (no), in other words 
by the discrete nature of the medium. 
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Figure 1: Representation of the selection rules i|18J) (V = 0.5) and graph of the 
solution lj20"|) . 



Inserting the dispersion relations (|16fl in the equation (|18|) . we obtain after 
some algebra the following solution 

2 sin(fci + K/2) =V, 2 sin(fc 2 + K/2) = -V. (20) 

These equations determine completely the wave nubers k 2 and K from the data 
of V (physics) and k\ (incident wave), as soon as one assumes that k\ an k 2 are 
close to, and less than tt. 

Remark : a particular case is obtained in the limit K — > for which k 2 — ► 
— k\ and 2 sin^) = V . This is the two-wave resonant interaction studied in 13 
for which a different multi-scale analysis has to be employed (resonant process 
occurs at larger scales). 

2.4 Three- wave resonant interaction 

Corresponding to the first selection rule l|18|l . the evolution equations for the 
envelopes that ensure vanishing of the resonant terms in the evolution l|13l) read 

[dt + 2smkid x ]A = —iBg, 

[d t -2sink 2 d x ]B = -iAg, (21) 

From now on we rename x\ by x and t\ by t. The same procedure applied to 
the second selection rule l|19[) eventually furnishes 

[dt + 2sin/ci9 a; ] J 4 = —iBg, 

[d t -2smk 2 d x ]B = -iAg (22) 

ft -v-=j»*. -«£>*. 

These are two standard 3-wave interaction nonlinear evolutions which are inte- 
grable systems on the infinite line x G M. when Cauchy initial data are prescribed 
in a space of functions vanishing (together with all their derivatives) at x — > ±oo 

na- 

Such an initial-boundary value problem is not the one we are interested in 
as indeed we study the scattering of an incident HF wave of envelope A(x, t) 
having a prescribed value for all time t at the origin of the medium, i.e. at 
x = 0. Moreover, as B(x,t) stands for the envelope of the reflected wave, its 
value will be prescribed to vanish at the output x = L for all time. Thus we 
cannot make use of the inverse scattering transform, unless first reformulated for 
a boundary-value problem on the finite interval, which is still an open problem. 

Note that the system ()22l) . by renaming g as ~g, maps to the first one except 
for the sign of the inhomogeneous term of the last equation. This change of sign 
is fundamental as it switches from instability to stability, as described below. 
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2.5 Stability properties 

The problem we consider is thus the scattering of an incident HF wave of 
envelope A(x,t), belonging to a carrier wave of frequency w\ % that generates 
backscattered wave of envelope B(x, t) and LF acoustic wave of envelope g(x, t) 
out of initial vacuum. Thus we perform a linear stability analysis of both sys- 
tems about the solution 

A(x,t)=A c , B = 0, 3 = 0. (23) 

with constant A c corresponding to a continuous wave (CW) irradiation. Math- 
ematically speaking this is an exact solution of both systems (|21[) and 121'l) and 
thus only an instability could produce an effective scattering. Let us seek now 
a solution as the perturbation 

A = A c + eae- wt 1 B = ebe~ iut , g = e q e~ lvt . (24) 

The system IJ21I) at order e then gives the linear system 



(25) 



possessing the 3 eigenvalues v = and v = ±i(\A c \/V)y/Q/2. The system is 
thus unstable (fl > 0), while the same analysis with system Ij22(l and g — eqe wt , 
furnishes the real eigenvalues and ±(\A c \/V)y / fl/2, and thus stability. 

Consequently the linear stability analysis predicts that the selection rule 
<|18|) will produce an effective scattering for an incident CW wave. Our purpose 
now is to demonstrate by numerical simulations that this instability is a soliton 
generator which induces an effective energy absorption of the incident radiation. 

3 Numerical simulations 

In order to understand the effect of the instability of the solution l|23(l in system 
H21|) . and compare it to system l|22|l . we perform here numerical simulations 
of those sytems under the following initial-boundary value data on the interval 
x€[0,L], 

A(0,t) = A c , B{L,t) = 0, g(x,0) = 0. (26) 

As an illustartion the figure [5] shows the energy density profile |A(a;,i)| 2 of the 
incident wave as a function of x at time t = 20 for A c = 0.7. The dashed line 
represents the acoustic phonon |g(cc,£)| 2 that has gained the energy lost by the 
incident optical phonon. 

This process corresponds to a brutal transfer of energy from the incident 
wave (which has already settled in the medium) to the medium, as shown on 
fig. |3 The normalized energy carried by the waves of envelopes A and g are 




Figure 2: Energy density at time t — 20 for optical phonon wave ^(a;, i)| 2 (full 
line) and |<7(a;,i)| 2 (cross). 
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Figure 3: Energy of the optical phonon wave I(t) for A c = 0.7 and the values 
(J2H) 



defined by 
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dx\A{x,t){' 



dx\g(x,t)\ 2 



(27) 
(28) 



These are the quantities evaluated on a numerical simulation in fig. [3] for the 
following parameters 



Wi = 3.94, fci = 2.89, 
uj 2 = 3.98, k 2 = 2.99, 
V = 0.1, ft = 0.04, if = 0.4. 



(29) 



Performing the same numerical simulations in the case of the system (|22H 
have never produced any nonlinear energy transfer, confirming the predictions 
of the linear stability analysis. 

4 Conclusion 

We have shown that a molecular chain allowing for wave coupling process, like 
the Davydov model, can present nonlinear energy absorption by resonant inter- 
action, coupling HF to LF waves, with selection rules (|18|l that are allowed only 
in discrete systems. 

The resulting governing equation, though being integrable (for a Cauchy 
initial value problem on the infinite line) leads to an instability when driven by 
boundary data on the finite interval. This instability is then the source of soliton 
formation in the medium, and energy transfers from the incident radiation to 
the medium excitation. 

This process is clearly illustrated by numerical simulations and opens the 
way to further studies in different contexts where the discreteness is known a 
priori to play a fundamental role. 
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